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Abstract 

We consider the scalar integral associated to the 3-loop sunrise graph with a massless line, two 
massive lines of equal mass M, a fourth line of mass equal to Mx, and the external invariant timelike 
and eqnal to the square of the fourth mass. We write the differential equation in x satished by the 
integral, expand it in the continuous dimension d around d = 4 and solve the system of the resulting 
chained differential equations in closed analytic form, expressing the solutions in terms of Harmonic 
Poly logarithms. As a byproduct, we give the limiting values of the coefficients of the (d — 4) expansion 
at a: = 1 and a: = 0. 


1 Introduction 


This paper is devoted to the analytic evaluation of the scalar integral associated to the 3-loop sunrise graph 
with a massless line, two massive lines of equal mass M, a fourth line of mass m, with m ^ M, and the 
external invariant timelike and equal to as depicted in Fig. |^. We use the same method as in jl], which 



Figure 1: The considered 3-loop sunrise graph. 


deals however with a 2-loop case. We start by writing the differential equation in a; = m/M satisfied by the 
considered integral (which turns out to be a third order differential equation), expand it in (d — 4), where 
d is the continuous dimension which regularizes the integral, and solve the resulting system of chained 
differential equations for the coefficients of the expansion in (d — 4) by repeated use of the variation of 
constants formula by Euler. The approach allows to obtain the coefficients of the expansion of the integral 
in (d — 4) up to virtually any order in (d — 4), expressing all of them in terms of Harmonic Polylogarithms 
(HPL’s 0) of argument x. At a: = 0,1 we recover the two quantities J 13 and Jn introduced in for the 
analytic evaluation of the {g — 2) of the electron at 3 loops in QED. J 13 , which corresponds to the vacuum 
graph of Fig.2, and Jn, which is the on-mass-shell value of the graph of Fig.l at a; = 1, where evaluated 
analytically in up to order (d — 4)^ and (d — 4)^ respectively. As a byproduct of the analytic result 
established for arbitrary a:, by taking its limiting values at a: = 0 and a: = 1 we obtain two more terms 
of the expansion in (d — 4) of J 13 (the analytic value of J 13 , exact in d, can be easily obtained by direct 
integration of the corresponding loop integral, so that this result can be seen as a check of the approach) 
and of Jn (a result not yet known in the literature). Those higher order terms in the (d — 4) expansion 
are expected to be needed in the analytic evaluation of 4-loop QED static quantities (such as {g — 2) of the 
electron, charge slope and renormalization counterterms). 


2 The differential equation 


Our aim is to evaluate the scalar integral 

F(d. M^ m*. = -m") = j 

In terms of the dimensionless variable x = m/M and putting P = Mp one can define the function 
as 


d‘^k2 d'^ks 

kim + M^){kl + M2)[(P -ki-k2- ^3)2 -b m2] ■ 


F{d,M'^,m'^,P^ = -m^) = C^{d) ^{d,x) , 


( 2 . 1 ) 
<i)(d, x) 
( 2 . 2 ) 


4 — d 

where C'(d) = ( 47 r)T“r (3 — d/2) is an overall normalization factor, with the limiting value 0(4) = 1 at 
d = 4. The previous equation can also be written as 


$(d, x) = 


d’^k2 


d'^k. 


C-Hd) [ _ 

(27r)3(^“2) J /c2(fc2 _|_ I)(fc2 _|_ _ ^2 - fc3)2 -|- a;2] 


{p^ = ■ 


(2.3) 
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Following iQl, for evaluating Eq. (|2.3| ) we will first establish the differential equation in x satisfied by 
$(d, cc). For writing the differential equation, we have to find the Master Integrals (Mi’s) of the graph of 
Fig.l . That is done by solving the Integration by Parts Identities (IBP-id’s) [ij, for the integrals associated 
with the graph of Fig. 1. We generated an initial system of 492 equations involving 449 unknown associated 
integrals. By using the symbolical solutions of the IBP-id’s it becomes a system of 205 equations in 76 
unknown integrals. By solving explicitly that last system, we were left with four Mi’s only, which we take to 
be the fully scalar integral of Eq. dP and three more integrals, having as integrand the same denominator 
as in Eq.(^) and numerators respectively equal to the scalar products (fci • ^ 2 ), (/c 2 • ks) and (P • ^ 2 ). Those 
four Mi’s satisfy a system of four I®‘ Order Differential Equations in the x variable. We rewrite them in 
terms of a single higher order equation for $(d, x). In so doing, it turns out that the fourth MI decouples, 
and d>((i, x) is found to satisfy the following 3’’'^ Order Differential Equation in x, exact in d, 

x^(l - x ^)^-3 <i)(d, x) -1- I [2 - (d - 4)] -I- [2 -I- 5(d - 4)] x^|x-^<i)(d, x) 

-b| [-6 - 8{d - 4) - 2{d - 4)2] -k [2 - 4(d - 4) - 6(d - 4)^] x^| x) 

_[8 + 14(d-4) + 6(d-4)2]x$(d,x) = ■ (2.4) 

Let us note that the above equation looks slightly simpler when written in terms of x^; but as x is the 
proper variable for expressing the final results, we use x since the very beginning. 


3 The T ^ 0 behaviour. 


Eor studying the x ^ 0 behaviour, we can try to solve Eq. (2.4) with the ansatz 


$(d, x) 




Eour possible values for the exponents are found: 


(3.1) 


ai = 0 , 

02 = -(d- 2) , 

03 = (2d - 5) , 

04 = (d-2), (3.2) 


where the last value is forced by the inhomogeneous term, while the first three behaviours are those of 
the solutions of the homogeneous part of the equation. By inspection of Eq.(^), we see that the aimed 
$(d, x) is finite in the x ^ 0 limit when d > 2. That is true even when the continuous parameter d is in 
the range 2 < d < 2.5; as in that range of d the exponents 02,03 give terms which diverge for x —> 0, the 
corresponding behaviours are ruled out. In the d —> 4 limit, that implies the absence of terms which for 
X —> 0 diverge as l/x^ (corresponding to 02 ) or behave as odd powers of x (corresponding to 03 ). That 
gives two condi tion s which can be used for determining two of the three arbitrary constants presen t in the 
solution of Eq.(^). Further, only even powers of x are found to appear inside the brackets of Eq. (^ t|) , so 
that the actual expansion of <i)(d, x) for x —> 0 has the form 


00 

$(d,x) = y~]a 2 fc(d)x 2 ^ 

k^O 

00 

+ ^ 52fc(d)x2'‘ . (3.3) 


2 





The equation determines the coefficients a 2 k{d) in terms of the first coefficient ao{d) (which must then fixed 
by some different method), while all the coefficients b^kid) are fully fixed by the inhomogeneous term of 
the equation. Explicitly, one finds 


02 (d) 


04 (d) 


(3d-8)(d-3) 
d (2d - 7) 


ao(d) , 


3 (3d-10)(3d-8)(d-3)(d-4) 

2 d (d + 2)(2d-7)(2d-9) ’ ’ 


do(d) 

b2{d) 


1 1 

“2 (d-2)2(d-3)(d-4)3 ’ 

1 1 
2 d (d-2)2((i-4)2((i-5) ’ 


(3.4) 


etc. 


4 The 


a; 


1 behaviour. 


A similar discussion applies to the a; —> 1 behaviour. When solving Eq. (^) with the ansatz 

$(d,x)=^(l-cr)d^ fEc«(d) (1-a:)") , 

one finds for the exponents Pi the values 


\n—0 


Pi 

= 0, 


P2 

= 1, 


Ps 

II 

to 

1 


Pi 

= 2, 

(4.1) 


the first three due t o th e homogeneous part of the equation, the fourth to the inhomogeneous term. 
By inspection of Eq.(^), one finds that the integral is regular at a: = 1, so that only the behaviour 
corresponding to Ps is ruled out. When expanded in powers of (d — 4) in the d ^ 4 limit, it would generate 
for X —> 1 terms in ln(l — x), which are therefore forbidden. That information, coming from the study of 
the behaviour at x = 1, provides with a further condition for the determination of the three integration 
constants. 


5 The expansion in (g? — 4 ). 

The 3-loop integral <i)(d,x), Eq.(|2.3D, develops a triple pole in (d — 4) when d —> 4, so that its Laurent’s 
expansion reads 

CX> 

$(x,d) = (d-4)" $(")(x) . (5.1) 

n——3 

, the coefficients (x) are found to satisfy the following system of chained 

differential equations 


By inserting it into Eq.(2.4 


~ d3 2/I _ 1 1 \ d /2 1 

_ dx3 yx"*"! — X 1 + xj dx2 \x2~*'l — x”*"l-l-x/dx yx"*"! — x 



$(")(a;) 
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1 + X J dx \x 

\n^^+^\x) , 


1 


1 


1 /2 


1 — a; 

2 

1 — X 
1 

1 — X 1 + X ) {n + Z)\ 


1 — X 1 + X 

1 1 


1 — X 1 + X 




(5.2) 


where = 0 when n < —3. Note that in the equation for the inhomogeneous term is just 

the proper term of the expansion in (d — 4) of the inhomogeneous term of Eq.( ^.4| ), which is known (as 
matter of fact, a rational factor), while in the equation for <i)(“^)(x), <i)(“^)(x) appears in the r.h.s. as a 
contribution to the inhomogeneous term and so on; given such a chained structure of the equations, it is 
natural to solve them recursively, starting from the equation for <i)(“^)(x) and going up to the successive 
coefficients of the expansion in (d — 4) of $(d, x). Summarizing, all the equations (E^) have the same form, 


_J. _L 

dx^ V X 1 — X 1 + X , / dx2 


- 2 


1 


1 


, - -4 - H- 

1 — X 1 + x/dx Vx 1 — X 


1 + X 


$(")(x) = 

K^^\x) , (5.3) 


where (x) stands for the inhomogeneous term, which can be taken as known when solving the equations 
starting from n = —3 and continuing in order of growing n. 


6 The solution. 

The homogenous equation associated to all the equations corresponding to different values of n is always 
the same, namely 


d3 


1 


■ 2 ( - H- — 

dx^ \x 1 — X l + xy dx^ 


-2(i 


1 


1 


1 — X 1 + x/dx 


4--411 


1 


1 


X 1 — X 1 + X 


(t){x) = 0 , 

( 6 . 1 ) 

so that all the Eq.s(^.3|) can be solved by Euler’s method of the variation of the constants. In this case the 
general solution can be written as 


4>("^(x) = +(/)i(x) 

+02 (a:) 

+03(a;) 


$(-) + y" ^y My)^3iy)^^My)My) 


in) r .7..0l(y)03(j/) -0'l(y)03(y) rAn)f^.^ 

2 J dy K [y) 




( 6 . 2 ) 


where the 0i(x), i = 1,2,3 are three linearly independent solutions of the homogeneous equation Eq.(|^), 
l+(x) their Wronskian and the = 1,2,3 the integration constants. 

Although algorithms for solving differential equations, even with simple rational coefficients, are not 
available in general, by trial and error it is not difficult to find (and in any case easy to verify) that Eq. (0) 
admits the three linearly independent solutions 

= ( 1 -^"), 

1 (1 — x^)(l — x^) 

“ 2 (^2 


01 (x) 
02 (x) 


4 


2(1 -x2)id(0;x) 









































+ Y^(1-3^") [i?( 0 ,-l;x)+i?( 0 ,l;x)] 
1 (x^ + 3)(3x^ + 1) 

“ ^ X ’ 


(6.3) 


where the functions i?(0; x), 77(—1; x), 77(1; x), i/(0, —l;x) etc. are harmonic polylogarithms i, i- The 
Wronskian of the three above functions (f>i{x) is 


W{x) = - 


(1 — x)^(l + x)^ 


(6.4) 


When inserting Eq.s(|6.3D, Eq.(^) and K^'^^{x), as defined by Eq.s(5.2,^), into Eq.(3.2), one sees that in 
analogy to previous work ||^ all the integrations can be carried out in terms of harmonic polylogarithms, 
up to any desired value of n, and one is left only with the problem of fixing, for any order n, the three 

(n) 

integration constants 4)) '. 

It is worth emphasizing that the quantitative determination of the integration constants can be obtained 
by simply imposing to the solution the qualitative behaviours for x —> 0 , x ^ 1 discussed in the previous 
sections Sec.(||^). More exactly, by imposing the right qualitative behaviour at order n, one obtains the 
explicite, quantitative value of the integration constants up to order n — 2 included (in other words: imposing 
the right qualitative behaviour to $(“^^(x) one fixes the first three constants and so on). We have 

evaluated the solutions up to n = 5 (the results involves harmonic polylogarithms of weight 8 ), and fixed 
the integration constants up to n = 3 included. The integration constants were found to be 


$1''^ = 0 


(for any k) 


= 

= 

= 

= 

.(3) 


7 

’ 

7 

% ’ 
3395 


= -- 


55296 128 

108691 97 


128 

$(-i) = 

2 3072 

‘T(3) ' 


(-a ^ (-2) ^ 

^ 256 ^ 

$i°^ = -^ + -C(3); 

^ 4608 24 ^^ ’ ’ 


37 

^ ’ 


1645 

36864 


960 


= — 


15149 


442368 


663552 

4181051 


768^^ ^ 5120' 


16 


C(5) ; 




256 *’*'^^ ’ 
1038085 


1225 


= -- 


7962624 9216 

56294669 19 


C(3) 


^,4 _ 117 . 5 . _ _ 
30720 256^^ ' 48 


5308416 


95 

1536 


C(3) 


10240 


63700992 12288 


2833 

18432 


C(3) 


81 


C(5). 


The explicit expressions of all the <i)*^"'^(x) which we have evaluated are too long to be listed here, and we 
give only the values up to n = 0 : 


$(“3)(x) 

$(“^^(x) 


24 12^ ’ 

+ — + —x^-^-x'^ — -x^ 77(0; x) ; 

96 32 96 8 ^ ^ ’ 


71 


7 X 


25 

384 
19 107 


35 


1 


1536 512 


1 


x 2 - 


1152 
559 
13824^ 


1 


1 


H—T ix 1 — -x'^ ) 77(0; x) — —x^ 77(0; x)77(0; x) 


16 


32 


12 


1 + 77(0;x) + ^xM 1 - ^x^ 77(0;x)77(0;x) 


1 ( 1 — x^)(l — x"*) 

m x2 


[77(-l; x)77(0; x) - 77(0; x)77(l; x) + 77(0,1; x) - 77(0, 


(6.5) 

( 6 . 6 ) 

(6.7) 


-l;x)] 
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( 6 . 8 ) 


H{Q-, x)i?(0; x)i7(0; x) + ^(1 - x^)[H{Q, 1; x)i?(0; x) - H{Q, -1; x)iJ(0; x)] 


+ 1 ( 1 - 1 =) 


H(0,0,-l;i)-H(0,0,l;i) + lc(3) 


Having the explicit form of the solution for any x it is almost trivial to obtain its value at x = 1 and x = 0. 
It is to be observed, in this respect, that the integration constants do not contribute to the 

X = 1 value up to order n, so that with the explicit value of the integration constants up to n = 3 included 
we can evaluate $(d, x = 1) in Fig.(2) up to order n = 5 included. By using the tables of we find the 



Figure 2: The considered 3-loop sunrise graph, in the case of equal masses, on mass-shell. 


following expansion in (d — 4) 


^(d, X = 1) ‘^='^ 

1 1 


7 


253 


2501 

'8((i-4)3 ^ ((i-4)2 “ 1152 (d-4) 13824 


+(d-4) - 


59437 


1 


+ (d-4f 
+ (d-4f 
+(d-4)4 


165888 18 


-h 


-H- 


2831381 
1990656 
117529021 
23887872 
4081770917 


2592 36 108^^ ^ 9 1080 

109403 


286654464 31104 


3115 2 , „ 21805,497 2 , 2 

~^ln2 - 


432 


C(3)-TrHn 2 

1296 ^ 108 


-ln '‘2 - 
9 

3053 
12960’ 
341 


16 

ya4 


— —ln‘^2-^04-^7r'^ln2 H-7r^ln^2 -|- —7r^C(3) H-ln®2 -|- - —(^(5) — 32a5 

54 9 180 9 36 ^ 15 12 ® 


+ (d-4)M- 


125873914573 


3439853568 


3386467 2 109403 2 

-TT — "" 


373248 


5184 


2 , „ 765821 


21805 


26789_ 

' 31104 ’ 

24211 

144 


3115, 3115 497 2 , 3„ 497 2..ox 

-ln^2 H- 04 -7r2ln^2-7rV(3) 

648 27 54 432 ^ 


1296 

3053 


7r^ln^2 


2160 


7r'^ln2 


C(5) 


568 


71, 5„ 

05-ln “2 -I- 

45 


41 


7r° -h 


10 


151^ 

1 ^’ 


7rHn^2 - + ^C(3)ln^2 


88 

y 


27 2723 

-7ry(3)ln2 -h 176C(5)ln2 - 176a5ln2 -f ln®2 -H-C(3)^ - 160 Og -H 176 

2 36 


+ 0{ {d-4f ) . 


(6.9) 


(For the definition of the mathematical constants, see for instance Table 1 of [Q). The terms up to (d —4)^ 
agree, within an obvious change of notation, with the formula for Jn given in 0; the last two terms, of 
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order (d — 4)'^, (d — 4)® are new in the literature, and are in full agreement with the numerical values for 
the same quantities given in 0. 

Somewhat surprisingly, the coefficients of the expansion of 4>(d, x = 0) in (d — 4) at order n are harder 
to obtain in this way, as they involve all the integration constants up to n included. On the other hand, 



Figure 3: The vacuum 3-loop graph. 


<i)(d, x) at X = 0 corresponds to the vacuum amplitude of Fig.(3), which can be easily obtained by direct 
integration of the corresponding Feynman graph amplitude, so that its value can be used in practice as 
a check on $(d,x = 0). Conversely, by using the explicit value of <i)(d,x = 0) obtained by the direct 
integration, one can fix all the integration constants by imposing the proper behaviour at x = 0 and the 
value at X = 0, namely ao(d) in the notation of Eq.( |3.3|) . 

For completeness, we list here the exact value of $(d, 0) = ao(d) 

d'^fci d‘^k2 d-^fca 


<i)(d, X = 0) = 


its expansion in (d — 4) is 

1 


c-Hd) 

(27r)3(^-2) 

1 


+ k2 + ks)^ 

r2(i (d 4)) F2(1 + i(d - 4)) F(-2 - |(d - 4)) 


32 (d - 4) r2(l - i(d - 4)) r{-2 - 2(d - 4)) r(2 + Ud- 4)) 


( 6 . 10 ) 


<i)(d, X = 0) = — 


1 


+ 


7 


1 


25 


1 


24 (d- 4)3 96 (d-4)2 

959 


384 (d - 4) 


1536 24*’*'^^ 


+ (d-4) 
+ (d-4)2 
+ (d-4)3 
+ (d-4)4 


6144 
10493 


7 1 

-C(3) 4- 

96^^ ^ 960 


-b 


24576 
85175 
98304 
610085 
393216 
1 


-7r'^C(3) 

960 ^ 


—C(3) - ^ 
384^^ ’ 3840 


C(3) + 
C(3)4 
C(7) 




1 


1536 

959 


3072 

1 


r4 _ ^^(5) +-, 

64^^ ^ 2688 


-C^(3) 

48^ ’ 


+ (d-4)5 -f 


4087919 

1572864 


C'(3) 
768^ ^ ’ 


6144 
83 

10493 
24576 
7 


12288 


-C(5) - ^ 

256^^ ’ 1536 


C(3)- 


959 

245760’ 


7r4^(3) _ —C(7) ■ 

3840 ’ 512^^ ’ 


15 

Itm 

13 


C(5)- 


TT® -b -C^(3) 

192^ ^ ’ 


25 


43008 


115200 


- Y^C(3)C(5) 


+ 0 { (d-4)6 ) 


( 6 . 11 ) 
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The above expansion corresponds and agrees (again up to trivial changes of notation) with the quantities 
Ji3 of 1] and /i4 of 0. 
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